AD  A118380 


TECHNICAL  REPORT  ARLCB-TR-82015 


NUMERICAL  SOLUTION  TO  AN  AUTOFRETTAGED  TUBE 
WITH  CONSTRAINING  WALLS  AND  END  CLOSURES 

Peter  C.  T.  Chen 


June  1982 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 

LARGE  CALIBER  WEAPON  SYSTEMS  LABORATORY 
BEN^T  WEAPONS  LABORATORY 
WATERVLIET,  N.  Y.  12109 


AMCMS  No.  611102H600011 


DA  Project  No.  1L161102AH60 
PRON  No.  1A2250041A1A 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


DTIO 

ELECTEI 


AUG  1  9  1982 


82  0?  1& 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other  author¬ 
ized  documents. 

The  use  of  trade  name(s)  and/or  manufacture(e)  does  not  consti¬ 
tute  an  official  indorsement  or  approval. 


DISPOSITION 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it 


to  the  originator 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whin  Data  Entered ) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

IS 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

NUMERICAL  SOLUTION  TO  AN  A'JTOFRETTAGED  TUBE 

WITH  CONSTRAINING  WALLS  AND  END  CLOSURES 

5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  authors; 

Peter  C.  T.  Chen 

6.  CONTRACT  OR  GRANT  NUMBERfa) 

>.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Armament  Research  &  Development  Command 
Benet  Weapons  Laboratory,  DRDAR-LCB-TL 

Watervlift.  NY  12189 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

AMCMS  No.  61 1 102H60001 1 

DA  Project  No.  1L161102AH60 
PRON  No.  1A225004 1A1A 

II.  CONTROLLING  OFFICE  NAME  AND  ADORESS 

US  Array  Armament  Research  &  Development  Command 
Large  Caliber  Weapon  Systems  Laboratory 

Dover.  NJ  07801  ,  ..  .  . 

12.  REPORT  OATE 

June  1982 

13.  NUMBER  OF  PAGES 

19 

1 4*  MONITORING  AGENCY  NAME  a  ADDRESS^/  different  Irom  Controlling  Olllce) 

IS  SECURITY  CLASS,  (ol  thla  report) 

UNCLASSIFIED 

ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  ( ot  thf  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  the  abelract  entered  In  Block  20,  II  dlttarant  from  Report) 

It.  SUPPLEMENTARY  NOTES 

Presented  at  1982  Army  Numerical  Analysis  and  Computers  Conference,  U.c  Army 
Engineers  Waterways  Experiment  Station,  Vicksburg,  MS,  3-4  February  1982. 

Published  in  proceedings  of  the  conference. 

IS.  KEY  WOROS  ( Continue  on  raratea  alda  II  naceeeery  and  Idantlly  by  block  number) 

Container  Autofrettage 

Elastic-Plastic 

Finite  Difference 

Gun  Tube 

20.  ABSTRACT  (C Mfeu  an  re  reran  aide  It  nmaaaaaty  mrd  Identity  Or  block  number) 

This  report  presents  a  numerical  study  of  a  container  autofrettage  process. 
This  process  uses  internal  hydraulic  pressure  to  expand  the  tube,  restraining 
containers  to  control  the  amount  of  tube  expansion  and  the  press  force  to 
hold  the  end  closures.  The  incremental  finite-difference  approach  developed 
recently  by  the  author  is  extended  to  obtain  numerical  results.  The  effects 
of  restraining  walls  and  the  press  force  on  the  displacements  and  stresses  are 
discussed. 


DO  I  JAM  71  1473  ED.TIOM  OF  »  NOV  66  IS  OBSOLETE  UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmt  Data  Bntarad) 


TABLE  OF  CONTENTS 


INTRODUCTION 

FINITE-DIFFERENCE  FORMULATION 

BOUNDARY  CONDITIONS  AND  INCREMENTAL  LOADING 

NUMERICAL  RESULTS 

REFERENCES 


LIST  OF  ILLUSTRATIONS 

1.  The  boundary  displacements  ua,  u^  as  functions  of  internal  10 

pressure  p  with  no  end  force  f. 

2.  The  radial  stress  distributions  during  loading  with  no  end  11 

force  f. 

3.  The  hoop  stress  distributions  during  loading  with  no  end  12 

force  f. 

4.  The  axial  stress  distributions  during  loading  with  f  -  0  13 

and  -0.6. 

5.  The  residual  stresses  due  to  complete  unloading  from  14 

different  stages  with  no  end  force. 

6.  The  boundary  displacements  ua,  u^  as  functions  of  internal  15 

pressure  p  with  f  ■  -0.6. 

7.  The  radial  stress  distributions  during  loading  with  f  -  -0.6.  16 

8.  The  hoop  stress  distributions  during  loading  with  f  m  -0.6.  12 

9.  The  residual  stresses  due  to  complete  unloading  from  different  18 

stages  with  f  **  -0.6. 


« 1 ; 


Accession  For 

NTIS  GRA&I 

□ 

DTIC  TAB 
Unannounced 

Justification _ 

By 

Distribution/ 
Avaiiability  Codes 
avail  and/or 
)ist  i  Special 


INTRODUCTION 


The  Importance  of  favorable  residual  stresses  In  an  autofrettaged  tube  la 
well  known. I  The  container  method  Is  one  of  the  autofrettage  processes 
currently  being  used  for  gun  tubes.  It  uses  Internal  hydraulic  pressure  to 
expand  the  tube.  Restraining  containers  or  dies  are  used  to  control  the 
amount  of  tube  expansion  by  means  of  a  small,  predetermined  clearance  between 
the  Inside  of  the  containers  and  the  outside  of  the  tube.  The  press  is  used 
to  simply  hold  the  end  closures  or  seals  In  the  ends  of  the  tube  and  to 
support  the  forces  of  the  Internal  pressure  on  the  closures. 

Many  methods  for  solving  the  partially  autofrettaged  problem  In  a  gun 
tube  have  been  reported. However,  the  effects  of  constraining  walls  and 
end  closures  on  the  residual  stresses  have  never  been  discussed.  This  report 
presents  a  numerical  study  of  the  container  autofrettage  process.  The 
finite-difference  approach  developed  recently  by  the  author^  Is  extended  to 


^Davidson,  T.  E.  and  Kendall,  D.  P. ,  "The  Design  of  Pressure  Vessels  For  Very 
High  Pressure  Operation,"  Watervliet  Arsenal  Report  WVT-6917.  Also  in 
Mechanical  Behavior  of  Materials  Under  Pressure  (edited  by  Pugh,  H.L.D.), 
Elsevier  Co.,  1970,  Chapter  2. 

^Hodge ,  P.  G.  and  White,  G.  N.  ,  "A  Quantitative  Comparison  of  Flow  and 
Deformation  Theories  of  Plasticity,"  J.  Appl.  Mech. ,  Vol.  17,  1950,  pp. 
180-184. 

^Chu,  S.  C.,  "A  More  Rational  Approach  to  the  Problem  of  an  Elastoplastlc 
Thick-Walled  Cylinder,"  J.  of  the  Franklin  Institute,  Vol.  294,  1972,  pp. 
57-65. 

^Chen,  P.  C.  T. ,  "The  Finite  Element  Analysis  of  Elastic-Plastic  Thick-Walled 
Tubes,"  Proc.  of  Army  Symposium  on  Solid  Mechanics,  1972,  The  Role  of 
Mechanics  in  Design-Ballistic  Problems,  pp.  243-253. 

^Elder,  A.  S.,  Tomkins,  R.  C. ,  and  Mann,  T.  L. ,  "Generalized  Plane  Strain  in 
an  Elastic,  Perfectly  Plastic  Cylinder,  With  Reference  to  the  Hydraulic 
Autofrettage  Process,"  Trans.  21st  Conference  of  Army  Mathematicians,  1975, 
pp.  623-659. 

^Chen,  P,  C.  T. ,  "Generalized  Plane-Strain  Problems  in  an  Elastic-Plastic 
Thick-Walled  Cylinder,"  Trans.  26th  Conference  of  Array  Mathematicians,  1980, 
pp.  265-275. 


1 


obtain  the  numerical  results.  The  material  Is  assumed  to  obey  the  von  Mlses' 
yield  criterion  and  the  Prandtl-Reuss  Incremental  stress-strain  relations. 

FINITE-DIFFERENCE  FORMULATION 

Consider  a  long,  open-end  thick-walled  cylinder  of  Inner  radius  a  and 
external  radius  b.  The  inside  surface  of  the  tube  Is  subjected  to  hydraulic 
pressure  p  and  and  end  force  (pira2)  is  applied  to  simply  hold  the  end  closures 
or  seals.  The  additional  force  f  on  the  end  closures  will  press  against  the 
tube.  The  amount  of  tube  expansion  Is  restricted  by  means  of  restraining 
containers  of  Inside  radius  C.  The  cross  section  of  the  tube  Is  divided  into 
n  rings  with  rj  -  a,  r2,...,  r^  »  p,...,  rn+i  -  b,  where  p  Is  the  radius  of 
the  elastic-plastic  Interface.  Since  the  material  behavior  is  nonlinear,  an 
Incremental  approach  is  used.  At  the  beginning  of  each  incremental  loading, 
the  distribution  of  displacements,  strains,  and  stresses  is  assumed  to  be 
known  and  we  want  to  determine  Au,  Aer,  Aeg,  Aez,  Aor,  Aog,  A oz  at  all  grid 
points.  According  to  the  Prandtl-Reuss  flow  theory,  the  Incremental  stresses 
are  related  to  the  Incremental  strains  by 

{Aoi}  «  [dij  ]  { Ae j  >  for  l,j  -  r,0,z  (1) 

and 

[dij]  -  2G[v/(l-2v)  +  6±j  -  Oi'oj’/B]  (2) 

where 

2  1  , 

2G  -  E/(  1+v)  ,  S  »  -  (1  +  -  H'/G)o2  ,  H’/E  -  u/(l-a)  , 

3  3 

°ra  *  (or+ag+oz)/3  »  °i'  ■  °i  “  %  , 

a  -  (1//2) [(ar-ag)2  +  (ao-az)2  +  (crz-ar) 2 ]  1/2  >  0q  (3) 
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E  Is  Young's  modulus,  v  Is  Poisson's  ratio,  is  the  Kronecker  delta,  aE  is 

the  slope  of  the  effective  stress-strain  curve,  and  oQ  is  the  yield  stress  in 

simple  tension  or  compression.  When  a  <  o0  or  do  <  0,  the  state  of  stress  is 
elastic  and  the  last  terra  in  Eq.  (2)  disappears.  Since  the  incremental 
stresses  are  related  to  the  Incremental  strains  by  Eq.  (1)  and  Au  -  rAeg, 

there  exists  only  three  unknowns  at  each  station  that  have  to  be  determined 

for  each  increment  of  loading.  Accounting  for  the  fact  that  the  axial  strain 
ez  is  independent  of  r,  the  unknown  variables  in  the  present  formulation  are 
(Ae0)i,  (Aer)i,  for  i  »  1,2,...  ,n,n+l ,  and  Aez. 

The  equation  of  equilibrium  and  the  equation  of  compatibility  are  valid 
for  both  the  elastic  and  the  plastic  regions  of  a  thick-walled  tube.  The 
finite-difference  forms  of  these  two  equations  at  i  ■  l,...,n  are  given  by^ 
[(ri+l”2ri)(di2)i  +  (-ri+i+ri)(d22)iK Ae0)i 
+  l(r1+1-2ri)(dH)i  +  (-r1+1+ri)  (d2i)i  ]  ( Aer)i 
+  ri(dl2)i+lUee)i+i  +  ri(dii)1+1(Aer)1+1 
+  [(r1+1-2ri)(di3)  +  (-ri+i+ri)(d23)i  +  ri(di3)i+i]Aez 

-  (ri+i-ri)(oe-or)i  -  ri[(or)i+1  -  (or)i]  (4) 

for  the  equation  of  equilibrium,  and 

(r1+1-2ri)(Aee)i  -  (ri+i-ri)(Aer)i  +  ri(Aee)1+1 

-  (ri+i-ri)(er-C0)i  -  ri[(ee)i+i  -  (ee)i]  (5) 

for  the  equation  of  compatibility. 


^Chen,  P.  C.  T. ,  "Generalized  Plane-Strain  Problems  in  an  Elastic-Plastic 
Thick-Walled  Cylinder,”  Trans.  26th  Conference  of  Army  Mathematicians,  1980, 
pp.  265-275. 
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BOUNDARY  CONDITIONS  AND  INCREMENTAL  LOADING 


The  three  boundary  conditions  for  the  problem  are 

(I)  (dl2)(Aee)i  +  (dn)i(Aer)i  +  (di3)iAez  -  -Ap  (6) 

(II)  (di2)n+l(Aee)n+l  +  <dl  l)n+l( Aer>n+1  +  (d13)n+lAez  “  0  (7a) 

before  contact  or  (Ae0)n+1  -  0  after  contact  ,  (7b) 

(III) 
n 

I  (ri+i-ri){ri[(d23)l(AE0)i  +  (di3)i(Aer)i)  +  ri+i ( (d23>i+l( AG0)i+l 
1-1 

n 

+  (d13)l+l(Aer)l+in  +  l  (ri+l“ri>lri(d33)i  *  ri+l(d33)l+l  1  Aez 

1-1 

-  Af  /  TT  (8) 

Now  we  can  form  a  system  of  2n+3  equations  for  solving  2n+3  unknowns,  (Ae0)i, 

(Aer)i,  at  1  -  1,2 . n,n+l  and  Aez.  Equations  (6),  (7),  and  (8)  are  taken 

as  the  first  and  the  last  two  equations,  respectively,  and  the  other  2n 

equations  are  set  up  at  1  -  1,2 . n  using  Eqs.  (4)  and  (5).  The  final 

system  Is  an  unsymmetrlc  matrix  of  arrow  type  with  the  nonzero  terms  appearing 
in  the  last  row  and  column  and  others  clustered  about  the  main  diagonal,  two 
below  and  one  above. 

In  order  to  Increase  the  efficiency  of  the  program,  an  adaptive  algorithm 
based  on  a  scaled  incremental-loading  approach  has  been  Implemented.  In  each 
step,  a  dummy  load-increment  such  as  Ap  is  applied  and  the  incremental  results 
Ao^  for  1  -  r,0,z  at  all  grids  are  determined.  For  all  grid  points  at  which 
0  "  I  l°il  I  <  o0»  we  compute  the  scaler  g's  by  the  formula 

g  -  l-  <r  +  [r*  +  4||4°1||2(o„2-||oill2>1/2>/||ao1||2  (9) 

where 

r  -  | | Oil |2  +  |  |Aoi| |2  -| | ot  +  Aoi| |2  (10) 
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(11) 


and  |  |c±| 1 ,  ||Aoi| |,  | |a±  +  A0l||  are  computed  by 

lloill2  -  -  l(or-oe)2  +  (O0-OZ)2  +  (cz-or)2) 

2 

Let  X  be  the  minimum  of  the  g's.  Then  X  is  the  load-increment  factor  just 
sufficient  to  yield  one  additional  point.  A  sequence  of  x(J)  can  be 
determined  for  all  steps  J  »  1,2,...  ,m  and  the  updated  results  are 

p(j)  -  p(J-D  +  x(j  )  Ap(J  ) 

o^(J)  »  +  X(J)Ao^-i)  ,  etc.  (12) 

NUMERICAL  RESULTS 

The  numerical  results  are  obtained  on  the  basis  of  the  following 
parameters:  a  «■  1.895",  b  ■  3.21",  c  ■  3.2275",  n  -  50v  E  ■  30xl06  psl,  v  - 
0.3,  oQ  ■  17xl04  psl  and  H*  -  0.  The  maximum  internal  pressure  applied  is 
13xl04  psi  (max.  p)  and  the  maximum  end  force  applied  against  the  tube  is  f  * 
-0.6pTi(b2-a2).  Introducing  the  dimensionless  quantity  f  -  f/[  n(b2-a2)p]  ,  we 
have  -0.6  <  f  <  0.  Since  the  end  force  required  to  simply  hold  the  end 
closures  or  scales  Is  -pxa2 ,  the  total  end  force  applied  on  the  end  closures 
Is  F  »  pxa2[-l  +  f(b2/a2-l)j.  In  order  to  discuss  the  effect  of  end  force  f, 
the  numerical  results  have  been  obtained  for  two  extreme  cases,  i.e.  f  ■  0  and 
-0.6. 

(a)  f  -  0.  In  this  case  the  total  end  force  applied  on  the  end  closures 
Is  F  >  -pxa2,  just  enough  to  support  the  forces  of  the  internal  pressure  on 
the  closures.  Therefore,  there  is  no  force  applied  at  the  end  of  the  tube. 

The  maximum  internal  pressure  (p  ■  13xl04  psi)  is  applied  Incrementally  in 
three  different  stages.  The  displacements  ua,  u^,  at  the  inside,  outside 
surface  as  functions  of  internal  pressure  p  are  shown  in  Figure  1.  In  stage 
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one,  the  elastic  solution  due  to  a  dummy  Internal  pressure  Is  applied  and  the 
scaled  factor  to  cause  initial  yielding  is  determined.  The  closed  form 
elastic  solution  together  with  the  von  Mises*  yield  criterion  may  be  used,  and 
the  pressure  factor  corresponding  to  Initial  yielding  is  p*/o0  ■»  0.36875.  In 
the  second  stage,  scaled  incremental-loading  approach  is  used  until  the 
maximum  allowable  outside  displacement  (c-b)  is  reached.  At  the  Instant  when 
the  contact  between  the  tube  and  container  first  occurs,  the  pressure  p/oQ  is 
0.57916  and  96  percent  of  the  tube  has  been  yielded.  In  the  third  stage, 
there  is  no  outside  displacement  and  internal  pressure  is  increased  in  20 
equal  steps  until  the  maximum  p/o0  -  0.76471  has  been  reached.  The  relation 
between  pressure  and  inside  displacement  is  almost  linear  in  this  stage  as 
shown  in  Figure  1.  The  results  of  the  displacements  at  the  end  of  three 
stages  are  represented  by  the  points  1,2,  and  3.  The  corresponding  results 
of  the  stress  distributions  for  or,  09,  and  oz  are  shown  in  Figures  2  through 
4,  respectively.  It  can  be  seen  that  the  stress  distributions  at  the  end  of 
three  loading  stages  are  quite  different. 

The  residual  stresses  after  unloading  completely  from  the  end  of  three 
loading  stages  have  also  been  obtained.  The  results  for  the  residual  hoop 
stresses  for  three  stages  and  the  residual  axial  stress  for  the  last  stage  are 
shown  in  Figure  5.  The  differences  in  residual  stresses  between  stage  two  and 
three  are  much  smaller  than  those  before  unloading.  Thus  a  further  increase 
in  internal  pressure  is  possible  in  the  presence  of  restraining  container,  but 
the  increased  pressure  makes  little  differences  in  the  residual  stresses.  The 
purpose  of  the  outside  container  is  to  prevent  occurrence  of  large 
displacements. 
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(b)  f  »  -0.6.  In  this  case  the  total  end  force  applied  on  the  end 
closures  is  F  ■  -pira2(0.4  +  0.6  b2/a2).  This  end  force  is  larger  than  that 
required  to  support  the  forces  of  the  internal  pressure  on  the  closures. 
Therefore,  the  end  force  applied  at  the  end  of  the  tube  is  f  ■  -0.6  pii(b:-a2). 
For  this  case  the  maximum  internal  pressure  is  applied  incrementally  in  four 
different  stages.  The  displacements  ua,  ufc,  at  the  inside,  outside  surface  as 
functions  of  internal  pressure  p  are  shown  in  Figure  6.  The  points  1  to  4 
represent  the  corresponding  results  at  the  end  of  each  loading  stage.  At  the 
end  of  the  first  stage,  initial  yielding  solution  has  been  obtained  and  the 
pressure  required  is  p  ■  0.34594  o0.  In  the  second  stage,  50  scaled 
incremental-loading  steps  are  applied  until  the  entire  tube  becomes  yielded. 

At  the  end  of  the  second  stage,  the  pressure  factor  p/o0  is  0.50194  and  the 
outside  displacement  is  still  smaller  than  the  clearance,  i.e.,  ujj  ■  0.86126  b 
0o/E  <  0.0175".  Since  the  material  is  assumed  to  be  ideally  plastic,  the  tube 
would  collapse  if  there  were  no  outside  restraining  containers.  A  very  small 
increase  in  internal  pressure,  say  Ap/a0  =  0.0001,  will  close  the  clearance 
between  the  tube  and  container.  The  instant  when  the  contact  first  occurs  is 
the  end  of  loading  stage  three.  After  the  contact  we  increase  the  internal 
pressure  in  29  equal  steps  until  the  maximum  pressure  lias  been  reached.  The 
relation  between  pressure  and  internal  displacement  is  approximately  linear  in 
this  stage  as  shown  in  Figure  6.  The  stress  distributions  for  or,  oq  at  the 
end  of  four  loading  stages  are  shown  in  Figures  7  and  8  respectively,  and  that 
for  oz  shown  in  Figure  4.  The  change  in  stresses  during  the  third  loading 
stage  is  too  small  to  be  shown  graphically  in  these  figures  but  the 
differences  in  displacements  are  large  as  shown  in  Figure  6. 
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The  residual  sCresses  due  to  complete  unloading  from  the  end  of  each 
loading  stage  have  also  been  obtained  and  some  of  the  results  are  shown  in 
Figure  9.  It  can  be  seen  that  the  differences  in  stresses  during  loading 
stages  three  and  four  are  quite  large  but  the  corresponding  residual  stresses 
are  very  close.  This  also  shows  that  the  effect  of  outside  containers  and  end 
forces  on  the  residual  stresses  is  small  but  their  effects  on  the  displacement 

*•'  I 

and  stresses  during  loading  are  large.  In  the  presence  of  the  press  force  on 
the  tube  end,  the  axial  stress  distributions  change  drastically  as  shown  in 
Figure  4  and  as  compared  with  the  case  of  no  end  force.  By  comparing  the 
results  for  the  residual  axial  stresses  as  shown  in  Figures  5  and  9,  we  can 
see  two  different  stress  patterns,  one  is  almost  the  reverse  of  the  other.  As 
a  result  of  extra  press  force  on  the  tube  end,  the  final  residual  stresses  can 
change  signs. 
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Figure  1.  The  boundary  displacements  ua,  uj3  as  functions  of 
internal  pressure  p  with  no  end  force  f. 
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Figure  3.  The  hoop  stress  distributions  during  loading 
with  no  end  force  f. 
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Figure  5.  The  residual  stresses  due  to  complete  unloading 
from  different  stages  with  no  end  force. 
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Figure  6.  The  boundary  displacements^  ua,  uj,  as  functions  of 
internal  pressure  p  with  f  =  -0.6. 
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Figure  9.  The  residual  stresses  due  to  complete  unloading 
from  different  stages  with  f  =  -0.6. 
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